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Abstract 

Based on the K-deformed functions (^-exponential and K-logarithm) and associated 
multiplication operation (^-product) introduced by Kaniadakis (Phys. Rev. E 66 
(2002) 056125), we present another one-parameter generalization of Gauss' law of 
error. The likelihood function in Gauss' law of error is generalized by means of the 
K-product. This K-generalized maximum likelihood principle leads to the so-called 
ft-Gaussian distributions. 
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1 Introduction 

Gaussian (or normal) distribution [1] is one of the most well-known and fun- 
damental distributions in many fields of science, e.g., an error distribution 
in measurement, a probability distribution of a fluctuating physical quantity 
in statistical mechanics, which accounts for the successes of Boltzmann-Gibbs 
exponential distributions [2]. For statistically independent iV events, it is well- 
known that the limiting stable distribution becomes a Gaussian one for suf- 
ficiently large N. This is the consequence of central limit theorem [3], which 
is the most famous theorem in mathematical probability theory. Historically 
however the Gauss derivation, as known as Gauss' law of error, played a very 
important role before establishing the central limit theorem. 
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Quite recently one of the authors et al. [4] have shown a one-real-parameter 
(g) generalization of Gauss' law of error, in which the g-generalized maximum 
likelihood principle leads to the so-called g-Gaussian distribution, by utilizing 
the Tsallis (/-deformed functions [5,6,7,8] and associated multiplication opera- 
tion, i.e., so called g-product [9,10]. Tsallis' entropy and g-deformed functions 
have been successfully applied in order to explain the ubiquitous existence of 
power-law behaviors in nature. Tsallis' entropy S q = (J dxf(x) q — 1)/(1 — q) 
is a one-real-parameter generalization of Boltzmann-Gibbs-Shannon entropy 
S BGS = — J dx f(x)\nf(x) and employed for generalizing the traditional 
Boltzmann-Gibbs statistical mechanics, thermodynamics, and Jaynes' infor- 
mation theory. 

The g-deformed exponential and logarithmic functions are defined by 



e Wq (x) = [l + (l-q)x} 1 + -\ (1) 

7-1-9 _ 1 

\n q (x) = ——, (2) 



respectively. Here q is a real parameter characterizing the deformation func- 
tions, and [x] + = max(x, 0). In the q — > 1 limit, the g-exponential and 
g-logarithmic functions reduce to the standard exponential and logarithmic 
functions, respectively. The following useful identity holds. 

exp g (a; + y) = exp q (x) ■ exp q ^ - - ^ _ . (3) 

As in the general usage for a Gaussian function, we here mean that any func- 
tion of the general form exp g (— ax 2 ) with a positive constant a as a g-Gaussian 
function. 

One of the most fundamental ingredients for proving the g-generalization of 
Gauss' law of error is the recently introduced g-product [9,10], which is defined 
by 

x l-q + y l-Q_l 
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On the other hand, Kaniadakis [11,12] has introduced another one-parameter 
deformed exponential function 

exp| K }(a;) = (vT+ k 2 x 2 + kx^J k , (5) 
and its inverse function 

/ \ / \ 

lll{K}(x) = Yk ' 

where n is a real parameter and takes a value in the range (—1, 1). Both k- 
deformed functions are symmetric under interchange of k «-> — k. In the k — > 
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limit, the ^-exponential and k- logarithmic functions reduce to the standard 
exponential and logarithmic functions, respectively. 

Note that the /t-logarithmic function can be expressed in terms of the two 
g-logarithmic functions with the different q indices as follows. 



1 1 



(7) 



However the ^-exponential function cannot be expressed in terms of the q- 
exponential functions. 

Note also that both q- and /t-exponentials shows asymptotic power-law be- 
haviors for large x: 



ex P g ( a; ) a; ~ 00 (( 1 -^) ;r ) 1 ^ > 



(9) 



whereas for a small value of x both functions behave as the standard expo- 
nential functions. Only in the intermediate region they are slightly different 
each other. Although the difference is not so large, this means that if an ex- 
perimental data is well fitted with ^-exponential function it should not be well 
fitted with g-exponential one, and vice versa. Until now in the literature there 
are only a few experimental evidences which are well fitted with ^-exponential 
probability distributions: the flux distribution of the cosmic rays extends over 
13 decades in energy [13]; the rain events in meteorology (the number density 
of rain events versus the event size) [13,14]; and the analysis of the fracture 
problem (the relation between the length of a transversal cut of the conducting 
thin ribbon and the electrical resistance) [15]. 

Kaniadakis [13] had already introduced a deformed algebra based on his k- 
deformed functions before establishing the concept of the g-product [9,10]. His 
^-deformed product [13] is defined by 

x <S> y = 

K 
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X 



+ 



y - y 
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y - y 
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which has the following properties 



associativity (x <g> y) ® z = x <g>(y <8> z) , 
unit element x®l — l®x — x, 

K K 

inverse element x®{l/x) = {l/x)®x = 1, 

K K 

and satisfying the following relations 



(11) 
(12) 

(13) 
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ln w (x®yj = ln w (:r) + \n {K} (y), 
exp {k} (x) (g) exp {k} (y) = exp {k} (x + y). 



(14) 
(15) 



In this way the ^-product also has suitable properties for generalizing Gauss' 
law of error. It is then natural to ask whether we can generalize Gauss' law 
of error based on the /t-product. This paper presents a positive answer to the 
above question. We show the ^-generalization of Gauss' law of error, in which 
the ^-generalized maximum likelihood principle leads to the /t-generalized 
Gaussian distribution. In the next section, we present the ^-generalized ver- 
sion of the maximum entropy principle in order to derive the /t-Gaussian 
distributions. In section 3, by utilizing the /t-product, the likelihood function 
in Gauss' law of error is generalized. We then prove the solution of the k- 
generalized maximum likelihood principle is the /t-Gaussian distribution. The 
last section is devoted to our concluding remarks. 



2 Maximum K-entropy principle 

Kaniadakis' K-entropy [11] can be defined [16] by 



S K = c K (l-J dxf(x) 1+K )+C- K (l- J dxfix) 1 '*), (16) 

With C K = ; r. (17) 

2k(1 + k) V ; 

S K reduces to Boltzmann-Gibbs-Shannon entropy S BGS — — J dxf(x)lnf(x) 

in the limit of k — > 0. 

Maximizing S K under the two constraints, 

/oo roo 
dxx 2 f{x)=a 2 , / dxf(x) = l, (18) 
-oo J —oo 



leads to 



Sf(x) 



(POO POO \ 

S K -P J^dxx 2 f(x)- 1 j_^dxf(x)J =0, (19) 



where (3 and 7 are Lagrange multipliers associated with the two constraints, 
respectively. It's solution is the so-called /t-Gaussian distribution, 

/(x)=exp w (-^ 2 - 7 ). (20) 

Note that 7 cannot be factored out from the argument of the /t-exponential 
function unless k — 0. This is one of the most different properties of the k- 
exponential function against the g-exponential function. In fact, applying the 
similar argument to the g-deformed functions, we can obtain the g-Gaussian 
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function f q (x) = exp q (—(3'x 2 — 7'), where two Lagrange multiplier f3' and 7' 
are introduced. Then, due to the property of Eq. (3), 7' can be factored out 

as 

f q (x) = ex Pg (-(3'x* - 7') = exp^-Y) ex P? . (21) 

In Fig. 1, /t-Gaussian functions are plotted for some different k values. 
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x 

Fig. 1. /-{-Gaussian functions Eq. (20) for some values of k with (5 = 1. For each 
curve 7 is determined so that the normalization is satisfied, k = corresponds to 
the standard Gaussian function. 



3 /t-generalization of Gauss' law of error 

Let us consider the same situation as conventional Gauss' law of error [1,4], 
i.e., we get n observed values: 

xi,x 2 , ■■■ ,x n e R (22) 

as the results of n measurements for certain observations. Each observed value 
Xi (i — 1, • • • , n) is a result of the measurement of identically distributed ran- 
dom variable Xi (i — 1, • • • , n). There exists a true value x satisfying the ad- 
ditive relation: 

Xi = x + e { {% = 1, • • • ,n) , (23) 
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where each of is an error in each observation of the true value x. Thus, for 
each Xi, there exists a random variable such that = x+Ei (i — 1, • • • , n). 
Every has the same probability density function / which is differentiate, 
because X ± , ■ ■ ■ , X n are identically distributed random variables (i.e., E ± , • • • ,E,, 
are also so). 

In order to prove the theorem for the ^-generalization of Gauss' law of error, 
we use the following lemma. Although the proof, which can be found in Ref. 
[4], is simple and compact, we here show it for the sake of being self-contained. 

Lemma Let be a continuous function from H. into itself and satisfying that 
T,i=i 4> (e») = for every n G N and e 1 , ■ ■ ■ , e n G R with £™ =1 = 0. Then 
there exists a£l such that (e) = ae. 

Proof. In the case n = 2, we can easily see that (— e) = —0 (e) for every 
e = ei,e2 G M. Moreover, in the case that n = 3 we have (ei + e<i) = 
0(ei) + (e2) for every ei,e2 G R. From this result and continuity of 0(e), 
it is easy to show that 0(e) must be a linear function of e, which prove the 
lemma. ■ 

Theorem For a given set of the data if the likelihood function 

L{ K } (6>) of a variable 9, which is defined by 

L {k} (9) = L {k} (x u x 2 , • • ■ , x n ; 9) 

= f(e-x 1 )®f(e-x 2 )®---®f(e-x n ), (24) 

takes the maximum at 

= r = -f>, (25) 

then the probability density function / must be a K-Gaussian distribution: 

/ (9 - Xi ) = exp {K} ( - a K (9 - x t ) 2 + C«) . (26) 

where a K is a /t-dependent positive constant, and C K is a /t-dependent normal- 
ization factor. 

Proof. Taking the /t-logarithm of the both side of the likelihood function 
L {k} (9) in Eq. (24) leads to 

ln w (L W (*) ) = E ln W (f(9 - X,)) , (27) 

where the property (14) is used. Differentiating the above formula Eq. (27) 
with respect to 9, we have 

^ ln w (L W (9) ) = g A l n{K} _ . (28) 
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n d 



ei = 8-Xi, (i = l,---,n), <f> K (e) = — ln w /(e) 



When 9 = 9* the likelihood function takes the maximum, so that 

0. (29) 

6=6* 

Let ti and </> K (e) be defined by 

(30) 

respectively. Then Eq. (29) can be rewritten to 

n 

i=i 

On the other hand, from Eq. (25) , we obtain 

n n 

E(**-*0 = E* = °- 



(31) 



(32) 



i=i 



i=l 



Now our problem is reduced to determining the function <f) K (e) simultaneously 
satisfying Eqs. (31) and (32). By using the lemma, <f> K (e) should be propor- 
tional to e, i.e., 

K (e) = -2a K e, (33) 

where a K is a /t-dependent positive constant, and the factor —2 is introduced 
for the sake of the simplicity of the final result. Integrating Eq. (33) w.r.t e 
gives 

ln {K} /(e) = -a K e 2 + C K , (34) 

where C K is a normalization factor which is /t-dependent through a K . Thus we 
obtain 

/(e) = exp {K} ^ - a K e 2 + . (35) 

A normalizable distribution requires that a K should be positive, and the nor- 
malization then determines the factor C K . 

Finally let us confirm that the extremum of the likelihood function at 9 = 9* 
is the maximum. By using Eq. (33), we obtain 



= -^{-2a K {9-x i )) =-2a K . (36) 
Since a K is positive, L{ K } (9*) is the maximum. ■ 
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4 Concluding remarks 

We have shown two different methods to derive the K-Gaussian distribution, 
i.e., both maximum entropy and maximum likelihood principles are general- 
ized by utilizing /t-deformed functions and associated multiplication operation 
(/t-product). Both ^-generalized methods reduce to the standard ones, respec- 
tively, in the limit of k — > 0. As the g-product is the most fundamental ingre- 
dient in the g-generalization of Gauss' law of error [4], the K-product enables 
us generalizing Gauss' law of error for the K-Gaussian. 

As seen in both q- and ^-generalizations of the maximum likelihood principle, 
we remark that for any set of one-parameter generalizations of exponential and 
logarithmic functions we can generalize the likelihood function by introducing 
a suitable multiplication operation or product. Let us here denote this product 
as ® a , where a stands for a real-parameter of the deformed functions. By using 
this product we can define the generalized likelihood function, which is similar 
to Eq. (24) but /t-products are replaced with the products (g> a . The associated 
Gaussian distribution can be obtained from this a-parameter generalization 
of Gauss' law of error. 

The difference among q-, k-, and other generalizations of Gauss' law of error 
apparently lies in the definition of each generalized likelihood function in terms 
of the generalized products. For the standard case, the standard products of 
each distribution f(6 — Xi) define the original Gauss likelifood function L(6), 
which is the k — > limit of Eq. (24). It is well-known that the standard 
products of / means all / are statistically independent. Unfortunately until 
now there is no consensus for interpretating any generalized likelihood function 
nor generalized product, although we feel that a certain kind of correlation is 
relevant [17]. 
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